Dynamics of resistive state in thin superconducting channels 
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We theoretically study how the dynamics of the resistive state in narrow superconducting channels 
shunted by an external resistor depends on channel's length L, the applied current j, and parameter u 
characterizing the penetration depth of the electric field in the nonequilibrium superconductors. We 
show that changing u dramatically affects both the behaviour of the current-voltage characteristics 
of the superconducting channels and the dynamics of their order parameter. Previously, it was 
demonstrated that when u is less than the critical value u cl , which does not depend on L, the 
phase slip centers appear simultaneously at different spots of the channel. Herewith, for u > u c i 
these centres arise consecutively at the same place. In our work we demonstrate that there is 
another critical value for u. Actually, if u does not exceed a certain value u C 2, which depends 
on L, the current-voltage characteristic exhibits the step-like behaviour. However, for u > u C 2 it 
becomes hysteretic. In this case, with increase of j the steady state, which corresponds to the time 
independent distribution of the order parameter along the channel, losses its stability at switching 
current value j aw , and time periodic oscillations of both the order parameter and electric filed occur 
in the channel. As j sweeps down, the periodic dynamics ceases at certain retrapping current value 
jr < jaw- Shunting the channel by a resistor increases the value of j r , while j 3W remains unchanged. 
Thus, for some high enough conductivity of the shunt j r and j aw eventually coincide, and the 
hysteretic loop disappears. We reveal dynamical regimes involved in the hysteresis, and discuss the 
bifurcation transitions between them. 

PACS numbers: 74.40.Gh, 74.81.-g, 74.78.Na, 74.20. De, 74.78.-w 
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INTRODUCTION 

Understanding of the mechanisms responsible for de- 
struction of the superconducting state is one of the most 
fundamental problems in physics of superconductors. It 
is known now that at the critical current a supercon- 
ducting sample generally does not simply enter the fully 
normal state. The building up of an intermediate state 
structure [l[ , flux flow [2( , and heating effects [3j, |4[ lead 
to complicated transition behaviour [5[. 

A number of very interesting phenomena have been 
observed for quasi-low-dimensional superconductors (see 
Ref. [6[ for a review) . The measurements of the current- 
voltage characteristics (CVCs) of such superconductors 
revealed a large transition region between the first oc- 
currence of non-zero voltage and the completely normal 
state 0, Q- Within this transition region the voltage 
increases in a scries of regular voltage jumps. The exper- 
imental and theoretical studies showed that every volt- 
age jump is associated with the appearance of a localized 
phase-slip center (PSC) [9|, [lCf in narrow superconduct- 
ing channels or a phase-slip line (PSL) [3, Is |ll| in wide 
superconducting films. Moreover, the electric properties 
of a superconductor between the superconducting and 
fully normal state is determined by a system of interact- 
ing PSCs [1(. Their dynamics is still poorly understood, 
and at the moment is of a great research interest |12H16|. 

Previously it was found out that the parameter u = 
T^p/rg, where r^, is the relaxation time of the amplitude 
of the order parameter (OP) and rg is the relaxation time 



of the phase of the OP, governs the dynamics of the PSCs 
in a narrow superconducting ring 12|, |13|. With this, the 
superconducting channel can demonstrate two dynami- 
cal regimes. If u < 1, then the phase of the OP relaxes 
faster than the amplitude, and PSCs occurs more or less 
simultaneously at different regions of the channel (see 
Fig. 3b of Ref. [H). In the case of u > 1 the PSCs ap- 
pear consecutively at the sam e sp ot of the channel (see 
Fig. 3a of Ref. [lj|). In Ref. [H| it was shown that the 
steps in the CVC can be associated with bifurcations of 
either steady or oscillatory dynamics of OP. These bifur- 
cations can substantially complicate dynamics of the OP 
and eventually lead to appearance of such sophisticated 
phenomena as multistability and chaos [14j. Thus, the 
appearance of multiple PSCs can be explained in terms 
of bifurcation phenomena. For example, Refs. 14H16| 
report that the period-doubling bifurcation can be in- 
volved in formation and evolution of multiple PSCs in 
the spatio-temporal dynamics of OP. 

Another interesting property of narrow superconduct- 
ing channels is that their CVCs are able to demonstrate 
quite wide hysteresis, which is found to arise from Joule 
heating and the strong temperature dependence of the 
resistance of the wire [3j, [l7|, [l8[ . Shunting of the chan- 
nel by a resistor narrows this hysteresis, and for small 
enough resistance the hysteresis disappears |18j . 

It was suggested that shunting is able to influence the 
phase-slip events, and hence can be used for control of 
dissipation and coherence in nanowires [18j. However, 
the dynamical mechanisms involved in the formation of 



the hysteresis are still not clear, and the ways of how 
shunting drives the phase-slip events are not fully appre- 
hended. In this paper we study the dynamics of the resis- 
tive state in a narrow superconducting channel shunted 
by an external resistor and explore the bifurcation mech- 
anisms involved in the hysteresis of CVCs and also re- 
sponsible for formation and evolution of PSCs. 

The article is organized as follows. In Sec. II we intro- 
duce the model, which we used to study the properties of 
the superconducting channels. In Sec. Ill we discuss the 
behaviour of the resistive state of the unshunted chan- 
nels for different values of the parameter u, and reveal 
the bifurcation mechanisms responsible for the destruc- 
tion of the periodic dynamics at the retrapping current. 
We study the superconducting channels shunted with an 
external resistance in Sec. IV. Finally, we summarize the 
results in the Sec. V. 



FORMULATION OF THE PROBLEM 

In the present paper we study the properties of a su- 
perconducting channel of the length L by simulation 
of the dimensionless time-dependent Ginzburg-Landau 
equations (TDGLEs) 0: 
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Here ip = pcxp(i6) is the dimensionless complex OP, 
where p is the amplitude and is the phase of the OP. 
The distance is measured in units of the coherence length 
£ and time is measured in units of the phase relaxation 
■■, A is the penetration depth, a n is the 
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normal state conductivity, and c is the speed of light. 
The electrostatic potential <fi is measured in the units 
of (/>o/27rcT0, where e is the electronic charge and h is 
the Planck constant, the dimensionless current density 
j is defined in the units of </>oc/87r 2 A 2 £. The penetra- 
tion depth of electric field Ie is an additional parameter, 
which characterizes the noncquilibrium superconductors. 
From Eqs. (U [2]) it follows that l E = £/u 1/2 . There- 
fore, the parameter u characterizes the penetration of 
the electric field in the nonequilibrium superconductors. 
Notably, the value of the parameter u introduced in the 
previous section depends on the material of superconduc- 
tor. 

Let us discuss the range of applicability of the TDGLEs 
in more details. One of the first attempts to analyse the 
simplified version of the TDGLEs similar to Eqs. (|l|2j) 
was discussed in Ref. [20(. In this work TDGLEs with 
u = 12 were used to describe the gapless case with large 
concentration of magnetic impurities t s T c << 1, where 



t s is the scattering time on magnet ic impurities, and T c is 
the transition temperature 2l|,|22j. Note that in the gap- 
less case the TDGLEs can be rigorously derived from the 
microscopic theory. In that respect they are exact and are 
valid in the broad range of temperatures. More realistic 
generalized TDGLEs were derived by Kramer and Watts- 
Tobin [23}. Their version of the TDGLEs with u = 5.79 
takes into account inelastic scattering by phonons and 
are valid for dirty superconductors with a finite energy 
gap TimT c << 1, where Ti m is impurity scattering time. 
In the the dirty limit the inelastic scattering time is large. 
It leads to the enhancement of the penetration depth of 
the electric field Ie- This process can be described by 
the generalized TDGLEs [24J assuming that u is small 
enough (u < 1). Note that the generalized TDGLEs 
have an additional time derivative term which can lead 
to different dynamics at j > j c [22j . However, as it is 
pointed out in Ref. [27| (Eq.(196) and Fig. 8, page 499), 



this term does not affect the dynamics of the solution 
at large currents, while in the vicinity of the critical cur- 
rent the influence of this term has quantitative character. 
Therefore, Eqs. (JTJ[2]) describe well the gapless supercon- 
ductors with large concentrations of magnetic impurities 
as well as superconductors with the gap near T c . On the 
other hand variation of the parameter u allows us to de- 
scribe qualitatively ordinary superconductors in the dirty 
limit in the large temperature range. Definitely TDGLEs 
can not be used to reveal the properties of the supercon- 
ducting channels at T — > 0. However, some qualitative 
predictions based on this model can be used to describe 
quantum phase slips in this limit 25|, |26| . 

In the present paper we model a thin superconduct- 
ing channel connected at the edges to the bulk su- 
perconductors. In this case, the boundary conditions 
are p(-L/2) = p(L/2) = 1 and d(f>(-L/2)/dx = 
d<fi(L/2)/dx = 0, where the phase 8 and modulus p of 
the OP satisfy the relation ip = pexp(i0). The absence 
of the electric field at the edges of the channel determines 
the gradient of phase d9(—L/2)/dx = dO(L/2)/dx = j. 



GAPLESS CASE VS. SUPERCONDUCTOR 
WITH A GAP 

In this section we study the superconducting channels 
of the lengths L = Lq, and L = 2Lq, where Lo/£ = 11-78. 
Here we use the Ginzburg-Landau depairing current for 
a channel of a finite length j c (L) as the reference point 
defining Jc {L )/j c = 1.016 and j c (2L )/j c = 1.002, 
where j c = 2/(3v3) is the Ginzburg-Landau depairing 
current for an infinite channel. 

The value of the parameter u plays the crucial role in 
the process of creation of the PSCs. Similarly to the re- 
sults of the previous research (l2j, |l3( , we found out that 
critical value of the parameter u, u c i = 1, corresponds 
to the change of the dynamics of the PSCs. Actually, for 



u > u c i the phase of the OP 9 relaxes faster then the 
amplitude p. It leads to the appearance of the PSCs at 
the same spot of the channel, one after another (Figs, 
la - lc). Our simulations suggest that such behavior of 
the PSCs remains qualitatively the same for any studied 
channel's length, and any studied value of the applied 
current. For u < u c ± the amplitude p relaxes faster then 
the phase 8. In this case the PSCs can emerge in dif- 
ferent regions of the channel (Fig. Id). Hence, one can 
distinguish three different ways for the appearance of the 
PSCs: 1. There is only one PSC and the amplitude of the 
OP reaches zero several times at the same location. This 
is the case of the consecutive PSCs; 2. There are more 
than one PSCs and the amplitude of the OP reaches zero 
at the same time at all locations; 3. There are more than 
one PSCs, but the amplitude of the OP reaches zero at 
a different time for each location. Hereafter we refer to 
the cases 2. and 3. as to the case of simultaneous PSCs. 
Besides u c \ , we reveal another critical value of the pa- 
rameter it, u C 2(L), which in contrast to u c \ depends on 
the length of the channel. The critical values, u c \ and 
u C 2, divide the parameter plane (L,u) into the regions 
corresponding to the different type of solutions of the 
TDGLEs (see Fig. 2). Dynamics of PSC for u < u^o, u c \ 
(domain SS in Fig. 2) was partly discussed in Ref. [14[. In 
this case the first PSC appears in the center of the chan- 
nel (see Fig. lc), and CVC in the vicinity of j = j c (L) can 
be approximated as V ex y/(j — j c (L))/j c (the dashed 



line in Fig. 3). Analysis in Ref. [14[ showed that the 
period of OP oscillations in the vicinity of j = j c (L) for 
u < u c i is described by formula 



T = 



ttV3(u + 2) 
2 3 / 2 u 



[ti - jc) /it 



1-1/2 
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Equation (|3]) (solid and dashed lines in Fig. 4b, which 
correspond to u = 0.5 and u = 0.8 respectively) demon- 
strates excellent agreement with the results of the nu- 
merical simulations (squares and triangles in Fig. 4b). 
The increase of the current leads to the period-doubling 
bifurcation, which implies a sudden increase of the pe- 
riod at (j — j c (L))/j c = 0.1 for L = Lq and u = 0.5, 
(J - ic{L))/Jc = 0.01 for L = 2L Q and u = 0.5, and 
U - Jc(L))/j c = 0.0046 for L = 2L Q and u = 0.8 (Fig. 
4), and corresponds to the change of the spatial position 
of the PSCs: two adjacent PSCs appear symmetrically 
with respect to the center of the wire (Fig. Id). Within 
the SS region in Fig. 2 it is also possible to find more 
complex solution of the TDGLEs, which are character- 
ized by larger number of the PSCs [14| . Remarkably that 
the appearance of every new PSC induces a jump in the 
CVC as shown in Fig. 3. 

The HC regime in Fig. 2 corresponds to the qualita- 
tively different behaviour of the CVCs. When j < j c {L), 
the steady-state solution is stable (see insets to Figs. 3a 
and 3b). However, the increase of j induces the volt- 
age drop along the superconducting channel. The criti- 




FIG. 1: The arrangement of the PSCs: a) in the vicinity of 
jsw for u > u c i; b) for j 3> jaw and u > u c i; c) after the 
period-doubling bifurcation; d) in the vicinity of j c (L) for 
u < Ud and in the vicinity of j r for u > w c i • 
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FIG. 2: Phase diagram showing regions of stability of different 
solutions of the TDGL corresponding to the step like CVC 
for the consecutive (SC) and simultaneous (SS) PSCs and 
the hysteretic one for the consecutive (HC) and simultaneous 
(HS) PSCs. Insets represent the CVCs corresponding to the 
different areas of the phase diagram. 



cal current j = j sw corresponding to this bifurcation is 
called the switching current [4| j c {L) = j sw . This switch- 
ing is associated with the appearance of a single PSC in 
the center of the channel (Fig. la). Further increase of 
j does not change the spacial position of the PSCs, but 
decreases the period of the TDGLEs solution. The OP 
is strongly suppressed in the small region in the center of 
the channel. The dimensionless amplitude of the OP in 
this area is p ~ 10~ 7 , while at the edges of the channel 
p = 1 (Fig. lb). The period of the solution tends to zero 
with increase of the current (see Fig. 4a and inset to Fig. 
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FIG. 3: The CVCs of the channels of the length a) L = L ; b) 
L — 2Lo for different values of the parameter it. The heavy 
solid lines represent the CVC of the norma l state. Th e dashed 
lines show the analytical formula (V oc ^/(j — j c (L)). Insets 
show a) hysteretic behaviour in the vicinity of j c (L); b) the 
CVC of the channels corresponding to the case u > 1. The 
arrows illustrate hysteretic behaviour of the CVC. 



FIG. 4: The period of the solution as a function of the cur- 
rent of the channels of the length a) L = Lq; b) L — 2Lo for 
different values of the parameter u. Heavy solid and dashed 
lines represent the analytical formula (Eq. 3) for u = 0.5 
and u = 0.8 respectively. Thin solid lines represent approx- 
imations by the logarithmic function (Eq. 4). The arrows 
illustrate the period-doubling bifurcation point. 



4b). 

If now one decreases the current, the superconducting 
state restores for j < j sw indicating of the presence of a 
hysteresis. The critical current, j = j r , below which the 
superconducting state reappears is known as the remap- 
ping current [Jj. The CVCs of the channels in the HC 
regime has a pronounced tail in the vicinity of j r (see 
inset of Fig. 3b). The finite effective resistance appears 
in the range of j < j r , which induces the change of the 
slope of the CVC. Similar effects were reported in [28( 
for capacitively shunted Josephson junctions. However, 



in contrast to the results of Ref. |28| , the system under 
consideration demonstrates a hysteretic type of CVCs. 
For j < j r the PSC appears in the center of the channel 
with significantly lower value of the voltage and higher 
value of the period comparing to the case of j > j r (Fig. 
lc). Such behaviour of the PSCs is similar to the case of 
u < u c i, when we observe one PSC in the middle of the 



channel in the vicinity of j c {L). 

The effective resistance of the channel grows with in- 
crease of the value of the parameter u (Fig. 3). However, 
for a large range of j it still remains much lower than 
the effective resistance of the normal state (solid line on 
Fig. 3). Nevertheless, with increase of the current the 
CVC of the superconducting channel in the resistive state 
approaches the CVC of the normal state. These results 
are in agreement with the experimental results [l|, |8|, |29| , 
where the CVCs of the superconducting channels at var- 
ious temperatures have been investigated. 

The hysteretic behaviour of the CVCs in the HC 
regime (Fig. 5a) can be illustrated by the transition be- 
tween different solutions of the TDGLEs. For j > j sw 
the TDGLEs has only one periodic solution, which is 
represented in the phase space by a stable limit cycle 
(Fig. 5b), and which corresponds to a single branch in 
CVC (see Fig. 5a). With decrease of j, at j = j sw a 
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FIG. 5: a) Hysteretic CVC corresponding to the HC domain 
on Fig. 2 obtained for parameter values L = Lq, u = 2; b)- 
f) Projections of the phase trajectories to [p(0),E(0)] plane 
illustrating the structural changes that occur along the point 
sequence a, b, c, d, and e in inset a) corresponding to the value 
of the current (j - j c (L))/j c = 0.05, 0, -0.016, -0.026, -0.03 
respectively. Full circle and cross represent the stable and 
unstable steady-state solution of the TDGLEs respectively. 



pair of steady states are born in the vicinity of the limit 
cycle as the result of saddle-node homoclinic bifurcation 
(Fig. 5c). With this, one of these steady states is sta- 
ble, whereas another one is unstable. The coexistence of 
two stable solutions (bistability), a steady state and the 
limit cycle, for j < j sw stipulates two branches in the 
hysteretic CVC illustrated in Fig. 2b and Fig. 5a. With 
further decrease of j the limit cycle approaches to the 
unstable steady state (Fig. 5d), and at j = j r they col- 
lide and d isap pear, suggesting a homoclinic bifurcation 
(Fig. 5e) [30|. Remarkably, that in the vicinity of this 
bifurcation the period of the limit cycle can be approxi- 
mated by the logarithmic law (thin solid lines on Fig. 4a 
and inset to Fig. 4b), similarly to the result obtained in 
Ref. 0: 



For smaller current j < j r the only limit solution of the 
system is the stable steady state, which corresponds to a 
single branch in the CVC in Fig. 2b and Fig. 5a. Thus, 



FIG. 6: The critical currents as a function of the parame- 
ter u for the channel of length L = Lq. Squares and circles 
represent retrapping and switching currents respectively. 



the destruction of the periodic solution at the retrap pin g 
current is associated with a homoclinic bifurcation [301 ]. 
Note that the similar way of disappearance of OP oscil- 
lations was described in [5j, [3l(. However, a particular 
instability leading to the destruction of the periodic so- 
lution was not revealed. 

The length of the superconducting channel L strongly 
affects formation of PSCs, and correspondingly CVCs. 
For example, in the case of the longer channels L > 2Lq 
the value of u c i drops below u c \ (Fig. 2). It leads to 
the unusual behaviour of the CVCs (see graph for u = 1 
in Fig. 3b) for u C 2 < u < u c \ (HS domain in Fig. 2) 
comparing with the CVCs corresponding to the other do- 
mains. Such unusual shape of CVCs can be understood 
in terms of evolution of PSC in the case of L = 2Lq, 
u = 1 (full circles in Fig. 3b). For these parameter 
values a single PSC (Fig. lc) appears in the center of 
the channel in the vicinity of j sw . Then, at the critical 
current corresponding to the period-doubling bifurcation, 
two PSCs appear symmetrically with respect to the cen- 
ter of the channel (Fig. Id). The CVC and the period of 
the solution behave hysteretically in the vicinity of the 
period-doubling bifurcation. Then, at higher currents, 
the system switches to the HC regime, where a single 
PSC (Fig. lb) appears in the center of the channel with 
the period lower than one of PSC in the vicinity of j c (L). 
Since the value of u c i in this case is close to u c i = 0.99, 
the value of the switching current (j sw = 1.002) is quite 
close to the value of the retrapping current (j r = 1.0019). 

In the case of the shorter channels L < 2L$ there is 
another characteristic range of u c \ < u < u C 2 (domain 
SC in Fig. 2). For these values of u the PSCs always 



appear in the center of the channel. Note that in this 
case the CVC of the channel also possesses the step-like 
behaviour (stars on Fig. 3a). However, the change of the 
slope of CVC is not accompanied by creation of a new 
PSC. The period of the solution also remains unchanged 
(stars on Fig. 4a) . Such behaviour of the CVC is similar 
to one in the HC regime, where a tail of CVC appears 
for j < j r . However, in the former case the change of 
the slope of the CVC is more smooth. The absence of 
the period doubling bifurcation for u r -\ < u < u C 2 is in 



agreement with the previous results 12j, where it was 
shown that for u > u c i the PSCs appear consecutively. 

The values of the switching j sw and retrapping j r cur- 
rents are equal in the SS and SC regime (see Fig. 6). As 
it has been shown above, exceeding the critical value u c2 
leads to the hysteretic behaviour of the CVC. It corre- 
sponds to the appearance of the difference between the 
critical currents (domain HC on Fig. 6). Increase of the 
length of the channel causes decrease of the critical value 
u C 2- When u C 2 becomes smaller than u c i, domain SC 
disappears. It corresponds to the appearance of the HS 
domain and shift of the SS domain to the region of small 
u, u < u c i- 

As it is clearly seen from Fig. 2, u C 2 sharply increases 
with decrease of the length of the channel L. It leads to 
the step like behaviour of the CVC for the very narrow 
channels. As it has been shown before, the CVC in the 
vicinity of the Ginzburg-Landau depairing current shows 
the square root behaviour in this region of the parame- 
ters 



14 



It is in the agreement with the earlier result 
of Aslamazov, and Larkin [32|, who have considered the 
CVC of the short (L < £) superconducting bridge. 



A SC CHANNEL WITH A SHUNT 

In this section we discuss the properties of a supercon- 
ducting channel shunted with external resistors of differ- 
ent conductivity. Shunting the superconducting channel 
leads to decrease of the current flowing through the su- 
perconducting cannel. Therefore total current is defined 
as: 
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(5) 



where a s h is the conductivity of the shunt, and elec- 
tric field in the shunt is defined as E = — (</>(L/2) — 
(j)(—L/2))/L. In our study we fixed the value of u = 5, 
which corresponds to a hysteretic CVC (see Fig. 3) for 
an unshunted channel (i.e. of a channel with a s h = 0). 
For u = 5 and a s h = a single PSC (Fig. 7a) appears 
for any values of j, however, increase of a s h can dramat- 
ically affect the periodicity of the PSCs. For example, 
for a s h = 3 a single PSC appears in the vicinity of the 
switching current during the up-sweep of j (Fig. 7a) . At 
the same time, down-sweep of j induces a second PSC 



-0.5 




FIG. 7: The arrangement of the PSCs in the case of a shunted 
channel with: a) one; b) two; c) three; d) four PSCs per 
period. 




(j-j c (L))/j c x10 2 

FIG. 8: The period of the solution as a function of the current 
of the channel of the length L — 2Lo for different values of 
conductivity. Inset represents the period of the solution for 
the value of conductivity a — 5. 



in the middle of the channel (Fig. 7b). Approaching the 
retrapping current j r increases the PSCs' number, and, 
finally, destroys the periodic solution. With this, every 
change of the periodicity of the PSC induces a disconti- 
nuity of the period of the solution of the TDGLEs (Fig. 
8). 

Growth of a s h is accompanied by the increase of j r , 
however, the value of j sw remains constant. Therefore, 
at certain critical value of a s h = <J C the value of j r and j sw 
coincides and the hysteresis disappears. Similar experi- 
mental results were obtained in Rcf. [18[. It was shown 
that shunting the superconducting channel with exter- 
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FIG. 9: The CVC of the channel of the length L = 2L for 
different values of conductivity of the shunt. The solid line 
shows the analytical formula (V oc y/(J — j c (L)). The arrows 
represent hysteretic behaviour of the CVC. 



rial resistors with different conductivities lead to different 
CVCs(seeFig. 2aofRef. 0). If the conductivity of the 
shunt exceeds a certain critical value, then the CVC has 
the step-like structure. In the opposite case the CVCs 
have hysteretic character, in agreement with our calcula- 
tions. It should be noted however that the independence 
of j sw upon cr s h in our calculations is inconsistent with 



the earlier experimental results 18[. This inconsistency 
may be associated with the neglecting of thermodynamic 
fluctuations in our model, whereas the values of j sw may 
be sensitive to the thermodynamic fluctuations. 

Above the critical value of the conductivity, the de- 
pendencies of the voltage V and the period of oscilla- 
tions of OP T on j become quite complicated (see Fig. 



9 and inset to Fig. 



For a s h > a c the dependencies 



V(j) and T(j) demonstrate irregularities for low values 
of j. These irregularities are associated with appearance 
of new patterns in PSCs formation characterizing by dif- 
ferent numbers of PSC. For high conductivity a s h ^> <x c , 
the irregularities in the dependencies V(j) and T(j) shift 
towards larger values of j. 

For example, for the channel of the length 2Lq the crit- 
ical value of a c = 3.562. Interesting that for a s h > 5 a 
final voltage along the channel correspond to appearance 
of two (not one!) PSCs in the center of the channel (Fig. 
7b), which is similar to the results published in Ref. [15[ 
Moreover, within a feasible range of j we were unable to 
find a regime with only one PSC in the channel. Increase 
of j upto 1.035 and 1.056 for a s h = 8 and 10 respectively 
induces a third PSC in the middle of the channel (Fig. 
7c), while further increases of j upto 1.1 and 1.14 for 
ash = 8 and 10 respectively adds one more PSC (Fig. 
7d). Remarkably, the PSCs for the above values of j 
are rapidly follow one after another with relatively long 



time interval, which separates next group of PSCs. The 
number of PSCs per period increases quickly with the in- 
crease of j, which finally leads to the irregular behaviour 
of the period of OP (Fig. 8). 



CONCLUSION 



In conclusion, we revealed the dynamical mechanisms 
involved in the formation of the hysteresis, and showed 
how shunting drives the phase-slip events. We also stud- 
ied the dynamics of the resistive state in a narrow super- 
conducting channel shunted by an external resistor, and 
explored the bifurcation mechanisms which are responsi- 
ble for formation and evolution of PSCs and which are 
also involved in the hysteresis of CVCs. 

We have investigated the current-voltage characteris- 
tics of the superconducting channels of a different length. 
We have found two critical values of the parameter u. 
They induce the change of the dynamics of the solution of 
the TDGLEs. One of them is independent on the length 
of the channel (w c i = !)• It governs the periodicity and 
spatial position of the PSCs. They can be found in two 
regimes: simultaneous and consecutive ones. The second 
critical value of the parameter u depends on the length 
of the channel (u = u C 2 (!<)). It changes the behaviour of 
the CVC from the step like to the hysteretic. As a result, 
the phase diagram showing the regions of stability of dif- 
ferent solutions of the time-dependent Ginzburg-Landau 
equations can be divided to four regions corresponding to 
the step like current-voltage characteristics and the hys- 
teretic one for the consecutive and simultaneous phase 
slip centers. The current- voltage characteristics of the 
short L <C £ superconducting channels possess the step 
like structure. The voltage in the vicinity of j c (L) is 
proportional to y/j — j c (L). We also have revealed that 
the destruction of the periodic solution at the retrapping 
current corresponds to the homoclinic bifurcation. The 
results of our calculations show the switching between 
two metastable states in the area of the high currents re- 
garding the switching current and u > u c i . We have also 
investigated the superconducting channels shunted with 
an external resistor. It has been shown that the width of 
the hysteresis loop decreases with increase of the value of 
conductivity and disappears at a certain a c . In the tran- 
sition from the hysteretic current- voltage characteristics 
to the step-like ones, some of them possess irregular be- 
haviour. Different arrangements of the phase slip centers 
have been found for both the unshunted channels as well 
as for shunted ones. It has been shown that the change 
of the slope of a current voltage characteristics leads to 
the appearance of a new phase slip center. 
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